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Legendre polynomials. This type of determinant has a 
probability background dealing in birth and death processes. 
The method of analysis used in this thesis is a new approach 
to solving this specific example. This investigation involves 
systems of differential equations and Prufer's analysis in the 
phase plane. The following are new results obtained in 
addition to solving the special case mentioned above. 
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possess Sturm-Liouville properties. As a special case, a 
different proof of Turan’s Inequality is provided. Finally, 
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I. INTRODUCTION 


A. THE PROBLEM 

The main purpose of this thesis is to investigate and 
establish the essential Sturm-Liouville properties for special 
eigenfunctions which are expressed in determinant form. In 
particular, a special case is presented where the elements of 


the determinant are Legendre polynomials such that 


(n>k) k=0,1,2,... 
(-1<x<1) n=1,2,3,... 


Pi(x) P,,, (x) 


Ga) 
Pane) Po ., (Xx) 


® x (X) ~ 








This type of determinant has a probability background 
dealing in birth and death processes and was studied by Dr. S. 
Karlin and Dr. J. L. MeGregor. In their researches, Karlin 
and McGregor encountered determinants of classical polynomials 
as examples. In order to complete their results, they needed 
to know such properties as: number of zeros, interlacing 
zeros, and completeness. Their obvious starting point was to 
construct an ordinary differential equation such that (1.1) 
was the solution. However, this differential equation was not 
ina form that suggested how to answer any of these questions. 
[Ref. 1] 

The method of analysis used in this thesis 18S a new 
approach to solving this specific example. The investigation 
that follows provides a constructive answer, not only to this 


Bt 


example, but, to a whole class of such problems. The 
following are new results obtained in addition to solving 
Karlin's problem. 

Special determinants of other hypergeometric functions 
also possess the essential Sturm-Liouville properties. 

As a special case of (1.1), a different proof of Turan's 
Inequality is provided in Chapter III. Turan's Inequality is 


of the form 


P(X), Poa) -1<x<1 


(1025 
Pasi (XX) Paye (X) 








n=0,1,2, * 8 8 


where the elements of the determinant are Legendre 
polynomials. 

Finally, several theorems are presented for 
Sturm-Liouville systems of differential equations with 
polynomial coefficients. In general, boundary value problems 
of systems are not well-posed. Until now, very little has 
been established concerning these systems with polynomial 


coefficients. 


B. STRATEGY FOR SOLVING THE PROBLEM 

The basic strategy developed here is to construct a second 
order differential equation such that (1.1) is a solution. 
Next, this differential equation is reduced to a special 
system of first order differential equations such that the 


2 


coefficients are polynomials. Then, the solution trajectories 
are analyzed in the corresponding phase plane by a method 
known as the "Prtifer analysis". From this analysis, several 
theorems of Sturm-Liouville type are developed for general 
systems of differential equations with polynomial 
coefficients. Finally, completeness and orthogonality are 
established with respect to a positive definite matrix weight 


function. 


Il. CONSTRUCTION OF A STURM-LIOUVILLE SYSTEM OF 
DIFFERENTIAL EQUATIONS 


A. DEVELOPMENT OF A NONSTANDARD STURM-LIOUVILLE DIFFERENTIAL 

EQUATION WITH SPECIAL DETERMINANT EIGENFUNCTIONS 

The objective of this section is to construct a second 
order linear differential equation, with appropriate boundary 
conditions, such that its solution is (1.1). For the general 
case, with arbitrary n and k, this method is somewhat lengthy 
and requires extensive manipulation of equations. However, 
considering a special case of k = 1, this method is not as 
complicated yet still illustrates all of the important steps. 
For completeness, the outline of the solution for the general 
case the same procedure and is provided in the Appendix. 

In the general case, the eigenfunction B(x) has the 


following form 


0.0 ee Seen, 9) (n>k) k=0,1,2,... 


P(X) Pry (X) 


®. (Xx) = (229 


(-1<x<1) n=1,2,3,... 








Setting k = 1 results in 


x % (3x*-1) 
@ 1 (x) = 
P, (xX) Pose) 


The strategy is to use recursion formulas and to develop 
three equations for $1(x), a GG) and oom as functions 
of P,(x) and Bax)". Then, from these equations, a system of 
equations in the matrix form AX = O is constructed where 
X # 0. Finally, the determinant of A is calculated and set 
equal to zero. This results ina differential equation having 
(2.1) as a solution. 

First, the equations for $.4(x), BCX), and an Tes) are 


developed. Substituting the recursion formula 


Paoy (2) =P, (xx) - AAI) Ph (x) 


moto (2.1) for Py.) (x) and Pi4) (x) results in 


P(x) xP, (x) - 2a) ph (x) 
® x (X) = 
Pa(x) xP, (x) - AXP) p(x) 


Since k = 1, P,(x) and Eex) are replaced by x and ll, 
respectively. In addition, by computing the determinant and 


collecting like terms, we obtain the following equation for 


B,C x)- 
(exe) _ x(1-x?) 
n,1"——5 — Py (x) a (2.2) 


Next, we differentiate both sides of (2.2). Thence, 
1, 1 (x) = AED ph (x) -xP, (x) - 22) Ph x) 


(2-3) 
és [PA] 





Then, we substitute the Legendre differential equation 


aie | (1-2) Pa (xe) |= (nea) P, (2) 


into (2.3) and collect like terms. Then, (2.3) becomes 


/ a (i= 1c) ' 
@ 5,1 (x) =(n-1) xP, (x) + “AS P(x). (2.4) 


Repeating the same procedure on (2.4) results in the 


following relation for B" Cx). 


®) 4 (x) = <B@et) + (n-1) P(x) + (n-1) xP4 (x) (2.5) 


Now, a system of equations in matrix form is constructed 
from equations (2.2), (2.4), and (2.5). The system takes the 


form AX = O, i.e., 


OF (x) =Btn")) + (n-2) (n-1) x a 


=10/ 2 (n-1) (1-x*) Hp, (x) | . 
AX =|, 3 (2) (n-1)x yrs F =O, (2.6) 
1-3 =x (lias 
Pa,1 0 ce a Pa 


Since the trivial solution can not be a solution of (2.6), A 
is not invertible. Therefore, the determinant of A must equal 
zero. Hence, by setting the determinant of A equal to zero we 
obtain a second order differential equation where $1 (%) is a 


solution. 


|A| = (1-x?) (14+3x?) @3. (2) - 6%(1-x?) 2 4(x) + 


(2.7) 
(n-1) (ames) x2+n-2h0,, (x0) = 0 
By setting A = n-1l, then (2.7) takes the form 
(1-x?) (14+3x?) @7 (0) - 6x(1-x?) Of (x) + ~— 


A\3 (143) xt+a-2 fo, , (x) = 0 


Note that, for each Legendre polynomial solution P,(x) to 
the Legendre differential equation, there is a corresponding 
logarithmic solution, Q,( x). This logarithmic solution obeys 
the same recursion formulas as P(x). Hence, by the same 


procedure, the second solution to (2.8) is 


P, (x) Sees 
0. (Xx) sOReeod 


Therefore, the general solution is 


x %%(3x*-1) x %%(3x*-1) 
C; sai Op. : (2.9) 
P(X) Peat) OX) = Ona, (x) 


Finally, translated boundary conditions on (2.9) are 
imposed so the second logarithmic solution is eliminated. 
Note that (2.8) is not in the standard Sturm-Liouville 
form. However, #, (x) still possesses all of the essential 
properties. 
B. REDUCTION TO A SYSTEM OF FIRST ORDER DIFFERENTIAL 
EQUATIONS 
The objective in this section is to reduce (2.8) to a 


system of differential equations of the special form 


(1-x?) &/ 


aD + (1+3x*) ¥ -1<x<1 
(1-x?) P/ y 


BO + Y? 


(2.10) 


where a, 8B, and y are polynomials in x, # = $, 0%), and ¥ is 
an arbitrary function in x. There are an infinite number of 
ways to transform (2.8) into a system of differential 


equations. However, by requiring polynomial coefficients we 


remove the singularities at t (osputs Hence, all solutions are 
well behaved on the interval -l<x<l. 

The strategy is to convert (2.10) into a second order 
differential equation and then equate the polynomial 
coefficients of 8”, 8’ and &. 


The system (2.10) has the following equivalent form 
(1-x?) (3x7+1) O”-| (aty+2x) a @’ - 


a! (3x241) ~asxrgGX tt) - Yeoh ® = 


First, note that the coefficients of 8% are the same. 


Next, equating the coefficients of 8’ and simplifying results 


in 
Y= -2x-a . (mia) 
Equating the coefficients of @ and substituting (2.11) for 
y leads to 
2 
A 3 (A+3) x2+A-2 | = -@!/(3x?+1) +a6x-8& ace 
: (2.12) 
, &(-2x-4a) (3x7+1) 


1-x? 


Dividing both sides by 3x*+t+l and simplifying yields 


A(3(A+3)x7+A-1]) _ _ ys , _@6x 
3x*+1 3x°4+1 
epi (2.13) 
(3x*+1) | a (-2x-a) 
ee oe 


Next, we consider components of (2.13) that have 3x*+l in 


the denominator. The idea is to select a so that 


A(3 (A+3) x7+A-1] -a6x = 0 mod(3x*+1) 


This condition on ais essential in removing the singularities 


at (oe For now, set a and 8 so that 


a@=px , &=q 


where p and q are constants. Substituting px, p, and q into 


(2.12) for a, a’ and B, respectively, and simplifying results 


in 
(-3A2-9A) x4+ (2A24+10A) x7+(A2-A) = 
(-3p?-9p-9b) x*+ (-p?+2p-6b) x*+ (-p-b) 
Next, we equate the coefficients. Thus, p = 2A and q = 
-A(A+tl1) which implies that a@ = 2Ax and 8B = -A(Atl1). 


Substituting back into (2.11), it follows that y = (A+1)(-2x). 


Therefore, the system (2.10) takes the form 
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2Ax® + (1+3x*) ¥ -1<x<1 
-A(At1)® + (A+1) (-2x) ¥ 


(1-x*) @/ 
(1-x?) ¥/ 


A simple substitution of (A+1)¥ = ¥ implies the following 


equivalent form 


(1-x?) ®’ = 2Ax@® + (A41) (14+3x7) B -1<x<1 
(1-x?) B/ = -AD + (A+1) (-2x) ¥ 


(2.14) 


Hence, (2.14) is a system of first order differential 


equations that is equivalent to (2.8). 


C. PROUFER ANALYSIS 

The Prufer analysis is used to study the oscillatory 
behavior of trajectory solutions in the corresponding phase 
plane [Ref. 2:p. 312]. It is a very powerful method that is 
extremely useful. The objective of this section is to use 
this method to establish the theorems listed below for general 
systems of differential equations, and then apply these 
theorems to (2.14). These theorems for systems are adapted 
from the analogous theorems for Sturm-Liouville differential 
equations [Ref. 2:pp. 312-353]. 
Theorem l. (STURM SEPARATION THEOREM) 

If (% ¥), and (# ?), are linearly independent solutions of 


the system 


ll 


and 


48(-Q) -[y-a@]?#>0 


where a, 8B, ¥, and Q are polynomials in x, then €, must vanish 
at one point between any two successive zeros of $,. In other 
words, the zeros of & and &, interlace. 
Theorem 2. (STURM COMPARISON THEOREM) 


Let (8 ¥) and (& ) be nontrivial solutions of the systems 


aw (8) = (9) (8) am o(8) - GB) (8) 200 


where a, &, ¥, Q. a, 8, ¥, and Q are polynomials in x, and 


48(-9) -[y-@]*>0 
0<B<8 4B(-0) -[y-a]*>0 
4(8-B) (- (0-0) ) -ly-Y- (a-@) ]2>0 


Then, & vanishes at least once between any two zeros of &, 
unless (& ) is a constant multiple of (8 #). 
Theorem 3. (OSCILLATION THEOREM ) 


Consider the associated Priifer differential equation 
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(1-x2) B= 1 6c0876+ [(A+1) y-Aa) sinOcos68- (A+1) Osin?6 


obtained from the system of differential equations 
g(1-x2)(®)) = (Ae (A#1)0) (®) sexe 
y/ AB (A+1)y] |B 


where a, 8, y, and Q are polynomials in x. If 


4£(-0) -[y-a]?>0 


and @(-1) = t, such that -x/2 < t < x/2 for each A, then the 
solution @ is a continuous and strictly decreasing function of 


A for fixed x on -l<x<l. Moreover 


lim @ (x) =-« aa lim 9 (x) =< 


= 00 —~—06o 


for -l<x<l. 
Theorem 4. Consider the system of differential equations 


(2.15) which satisfies the Oscillation Theorem. If 


@(-1) <t, 5 <n for al]-<i<@ 
6(1)=o-nt1, = «oc n=0,41,242,... 


then (2.15) has an infinite sequence of real eigenvalues 
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soja Sg SA 


with 


lim A,==° ang 11M An=-° 


n-00 Nn-+-© 


The eigenfunction &, belonging to the eigenvalue A, has exactly 
n zeros in the interval -1l<x<l, and is uniquely determined up 
to a constant factor. 
1. Phase Plane 
The objective of this section is to convert to polar 


coordinates the following general system of differential 


equations 
-x*7)@' -= a@ + OF 
Cae . ad -1<x<1 (2.16) 
(1-x?)¥/ = £8 + yy 
where a, 8, ¥, and Q are polynomials in x. By introducing 


polar coordinates 


®=r(x)cos6(x), P=r(x) sin (x) 


where ® is positive in the counterclockwise direction, the new 
dependent variables become r and @. 
Expressing r and ®@ explicitly as functions of & and ¥ 


leads to 


14 


lie) 


@-arctan( at (2.17) 


y*=O2+ 72 . (2.18) 


Then, by differentiating both sides of (2.17) and 
(2.18), an equivalent system of differential equations for 


(2.16) is presented as functions of r and @. Differentiating 


both sides of (2.17) yields 


sec? 6 = - ¥e-0'¥ 


92 (2.19) 
Since 
: : ¥\* o2+9? 
sec’@ = 1+tan76 = 1+($ = ae 
then (2.19) becomes 
@® — 
& . at Fi ; (2.20) 


By substituting from the general system (2.16), 


equation (2.20) takes the form 


@ _ [£b+y¥) O-[ad+0F] ¥ 
(1-x") ae ¥ res oF 
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Substituting in polar coordinates, the differential 


equation for ® becomes 


(1-22) B=gc0s70+ (y-a) sin6cos8-Osin*6 (2.2%) 


Finally, we differentiate both sides of (2.18) and 
follow similar steps. The differential equation for r then 


becomes 


(1-x*) rer [acos*@+ (O+8) sin@cos8+ysin76] (2.22) 


Therefore, (2.21) and (2.22) represent an equivalent 
system to (2.16). 

The advantage of this system over (2.16) is that 
(2.21) isa first order differential equation of only ®. This 
allows the study of the oscillatory behavior of the solution 
trajectories independent of r. To find r for any solution 86, 
one solves (2.22) directly. 

We also note that ris always positive for nontrivial 
solutions. This can be seen by examining the solution of 
(2.22) [ReEE. 2:p 235277 

Next, we recall the following result from elementary 


algebra. A function 


F(x,y) = a x* + bxy +c y’* 


is negative definite iff and only if 
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as0O and 4ac - b*>0 


and positive definite if and only if 


a0 and 4ac - b*>0 , 


Next, we apply this result to (2.21). If we can show 


that 


48(-Q) -[y-a]7>0 


then equation (2.21) is either negative definite or positive 
definite, depending on the sign of 8. If equation (2.21) is 
positive definite, then the change in the angle ® is always 
greater than zero. This means that the angle @ is always 
increasing. Since our polar coordinate system is set up so 
that the positive direction of rotation is counterclockwise, 
then the solution trajectory of @ always rotates 
counterclockwise. If equation (2.21) is negative definite, 
then the solution trajectory of @ always rotates clockwise. 
In addition, since the change in the angle @ is either 
positive definite or negative definite, then the direction of 
the solution trajectory @ can never change. Therefore, we can 
conclude that once a solution trajectory crosses an axis it 
can not cross back. 

Also, since & = rcos@, then the zeros of $ occur when 
@ = x/2+nx where n = 0, +1, #2, 
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Consider the system of differential equations (2.14). 
The associated Pruifer equation is negative definite for A>O, 
and positive definite for A<O. However, for (2.14), it is 
easy to determine the direction of rotation across the &%-axis 
and ¥-axis without the Priifer equation. First, we examine the 
case when A>0O by directly substituting into (2.14). 

By setting @ = 0 we can determine the direction of 
rotation along the ¥-axis. Also, by setting ¥ = O we can 
determine the direction of rotation along the $-axis. Hence, 
if we substitute these values back into (2.14), then the 


corresponding phase plane for A>0O takes the form in Figure 1. 


Figure l. Phase Plane Analysis Along Y-axis and &%-axis 
for A>0O 
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Clearly, the solution trajectories can cross the @- 
axis and ¥-axis only in a clockwise direction. In the case 
when A<0O, the trajectories rotate in the counterclockwise 
direction. 

2. Sturm Separation Theorem 
Suppose (% ¥), and (% ¥), are two linearly independent 


solutions of the system 


(1-x?) ($"] (s 4 (=) -1¢x<1 (2.23) 


and 


48(-Q) - [y-a@]2#>0 


where a, 8B, y, and Q are polynomials in x. 

There are two cases to consider, one being when B<0 
and the other when 8>0. First, we consider the case when B<O. 

We examine the area "A" of the triangle created by 
joining the points (8 ¥),, (& Y),, and the origin in the (#,f) 
plane, Figure 2. 

By inspection, the area "A™ is equal to zero only when 
(@ ¥), and (8 ¥), are on the same straight line through the 
origin. This can occur with (# ¥), and (# ¥), both in the same 
quadrant, or in opposite quadrants. However, this can never 
occur Since (& ¥)) and (& YT), are linearly independent, so "A" 


has positive area. 
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Figure 2. Area Principle Analysis 


Recall that two nontrivial solutions of (2.23) are 
linearly dependent if and only if they are located on the same 
straight line through the origin in the (8,¥%) plane. Hence, 
the angle, o, between the linearly independent solutions, 

(s Y), and (@ Y),, is always O0<o<x. In other words, (8 ¥), and 
(& Y), are always in the same rotating half plane and never 
cross each other. 

Analytically, we can show that the area "A" is never 
zero by the following method. Using the area principle, we 


calculate the area "A" using the following formula 
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= o(@,¥,-@,%,)  . (2.24) 








The change in the area is 


= = ¥2(0)¥, + ®,%; = ©:f, = .¥5] 


Substituting in §8',, &',, ¥',, and ¥', from (2.16) and 


Simplifying gives 


dA 
Fe that) (OF, - OF) 


Substituting into this equation from (2.24) results in 


Aa(aty)a 


Finally, by solving for "A", the following formula is derived 


to determine the area "A" at any point -1l<x<l1; 


A-k exp| [iawn ae} 


where k = A(-1) is the initial area. 
Since a and y¥ are polynomials and -l<x<l, the area is 


finite. Moreover, A(-1) can not equal zero since the two 
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solutions are linearly independent. Hence, the area "A" can 
never be zero. 


Next, we recall that B<0O and 


48(-Q)-[y-a]4>0 


Therefore, d@/dx is negative definite which implies that all 
solutions rotate in a clockwise direction. Also, once a 
solution crosses the Y-axis, it can never cross’. back. 
Finally, since the area can not equal zero, the angle @ 
between (4 ¥)) and (84 Y), is always satisfies 0<o<x. 
Therefore, as the solutions rotate clockwise in the (#,#) 
plane, (& Y), and (8 Y), always cross the Y-axis in the same 
order. Therefore, the zeros of # and & interlace. The 
analysis for the case when B>0 proceeds in exactly the same 
way except that the solutions rotate counterclockwise. q.e.d. 
3. Sturm Comparison Theorem 
Let (& ¥) and (&€ n) be nontrivial solutions of the 


systems 


respectively, where a, B, ¥, Q, @, B, ¥, and Q are polynomials 


in x, and 
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46(-09) -[y-a]*>0 
o<B<B& 4E(-0) - [y-a@]?>0 
4 (8-8) (- (0-0) ) -[y-y-(a-a@) ]7>0 


Again, there are two main cases to consider. One case is when 
B-B > 0 and the other when B-B < 0. The second case is the 
same as the first except that the solutions rotate in the 
counterclockwise direction. Hence, it suffices to consider 
only the case when B-8 > O and that all solutions rotate 
clockwise. The objective is to analyze the two solutions in 
the (%,¥) plane through the use of the Priifer equation. 

The associated Prutifer equations for (2.25) are listed 


below, respectively. 


(1-x2) St =6c08%, + (y-a) sin®,cos, - Qsin’@, (2.26) 


(1-x?) = = Bcos70,+ (y-a@) sin®@,cos@, - Osin’@, (2.27) 


By subtracting (2.27) from (2.26), the resulting 
equation indicates the rate of change of the angle formed 
between (8 ¥) and (&€ nn). In addition, by superimposing one 
solution on top of the other, with 6, = 6, (Figure 3), then the 


equation becomes 
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(1-x*) 2 (8, -8,) = (8-B) cos76, + [y-y- (a-a@) ] sin8,cos®, 


-(Q-0)sin70,  . 


The resulting equation is a Priifer equation. This 
equation facilitates the analysis of the solution vectors 
(@ ¥) and (& n) by examining the behavior of the two angles, 
®, and @,. 


Since B-B > 0, and 


4 (B-8) (- (0-0) ) -[y-y-(a-a) ]7>0 


then, d/dx(@, -6,) > O. This means the rate of change of the 
angle between (% ¥) and (€ n) is always increasing, and (&€ ) 
"runs away'’ from (8 ¥), Figure 4. 

Note that (& ¥) and (& n) are solutions of two 
different systems of differential equations. The two solution 
may be linearly independent or dependent. Hence, one solution 


could "lap" the other solution one or more times. 
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Figure 3. Superimposition of Solutions 


~ 


(¥) 





Figure 4. One Solution “Runs Away" from the Other 
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Recall that the zeros of & and & occur when (8 ¥) and 
(§ ) respectively cross the Y-axis. Therefore, as the 
solutions continue to rotate, between any two zeros of # there 
is at least one zero of &. Note that, if (% ¥) = c(& n), for 
all -l<x<l, then the zeros occur simultaneously. Hence, & 
vanishes at least once between any two zeros of &, unless 
(— n) is a constant multiple of (% ¥). q.e.d. 

4. Oscillation Theorem 
We consider the system of differential equations of 


the form 


®/ A A+1 & (2.28) 
(1-x*) ($/) = (is (31) 9) ($) -1<x<1 


where a, B, ¥, and Q are polynomials in x, and 


48(-9)-[y-a]*>0 


The system (2.28) is basically the same as (2.23) 
except now we have introduced an eigenvalue, A. The objective 
now is to vary A thereby to study the zeros of &. The Prtifer 


equation for (2.28) is 


ae 
(1-x?) —+1 =) &cos?6,+ [ (A+1) y-Aa] sin8,cos8 
dx 1 Y i i (2.29) 


~(4+1)Qsin?6, . 
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As before, the zeros occur only on the ¥-axis which is 


the same as saying the zeros occur at @ = x/2+nx for n= 0, 
goles t2,..%- Next, we impose endpoint conditions on @ such 
that 
6(-1) =t, 5 KKK S for al]l-»<)<o 
(2.30) 
@(1) =o-nx, SKS n=0,41,42,... 


so that the angle t is as close to x/2 as possible. Note that 
the endpoint conditions are computed from tan@ = ¥/8. 


Now, we consider a second system 
2 E/\ _ (pa (pt+1)o) (Ee) _ 
(1 x)(¥)) (ue (h3) 9) (5) 1<x<1 


where p is an eignevalue. The Pruifer equation for this system 


is 


ca.) i =p6cos*68,+[ (+1) y-pa} sin8,cos6, 
(2.31) 


- (p+1) Osin’6. 
Suppose that x is fixed and the differential equation (2.31) 
satisfies the following endpoint condition 


6(-1) =t, —<«Ks for <p <eo 
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Following the same basic proof as in the Sturm Comparison 
Theorem, we subtract (2.29) from (2.31), and equate 8, to 6, 
at @)(-1) = @,(-1) = t. Hence, 


(1-2) 2 (6,-8,) = 
dx “* (2.32) 


(p-A) [Bcos76, + (y-a) sin8,cos8, -Qsin*6,] : 


Since 


4£(-Q) -[y-a]7>0 


then the discriminant of (2.32) is 


(p-A) *[48(-Q) - (y-a) 27) >0 ; 


Therefore, d/dx(@, -@,) is either positive definite or 
negative definite. Accordingly, there are several cases to 
consider. However, by examining one case carefully, all other 
cases follow directly, with only the direction of rotation 
varying. Thus, it suffices to show only one case. 

Let x be fixed. Suppose that 0<p<A which implies that 
p-A< 0. Now, consider the case when B<0O. 

Since B<0O and p-A<O, then d/dx(®, -@,) > 0. Hence, by 
considering the conditions that O<p and O<A, it is implied 
that all solutions rotate in the clockwise direction, and 
($ ¥) “runs away" from (& ) as A goes to positive infinity, 


Figure 5. 
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Figure 5. (# ¥) “Runs Away" from (& ), A>0 


If A<O, then p-A > O which implies that 
d/dx(@,-8, ) <0. Therefore, when A<0O, all solutions rotate in 
the counterclockwise direction and (% ¥) "runs away" from (& 
n) as A goes to negative infinity, Figure 6. 

In summary, as A goes from -® to », the angle 8, goes 
from © to -® in a strictly decreasing manner. Hence, for 
fixed x, @ is a strictly decreasing function of the variable 


A. In addition, 


Lan Q=-c and lim OQ =00 


A--2 


q.e.d. 
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Figure 6. (# ¥) "Runs Away" from (& n), A<0O 


5. Sequence of Eigenvalues and Interlacing Zeros 
Consider the system of differential equations (2.15) 
which satisfies the Oscillation Theorem. Now we continue with 
our analysis and maintain the conditions of the previous 


proof. Namely, let O<p<A and 


@(-1) =<, Bees. for all -~<j<@ 


6(1)=<o-nnz, > <<, n=0,41,22,... 


where t is as close to x/2 as possible. 
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Now, we let wp be fixed and 8<0. Note that a value for 
A which satisfies both endpoint conditions is an eigenvalue 


for the system of differential equations (2.28). Let @,(-1) 


= T. Then, as A increases, the solution trajectory of 8, 
rotates in a clockwise direction. Since d@/dx < 0, the 
trajectory of ®, eventually crosses @(1) = o-nx (n = 0). 


Hence, the first time that the trajectory of e) crosses the 
line @,(1) = o-nx is when n = 0. We call this eigenvalue Ag: 
Therefore, when 4 = A), the solution trajectory of ®, does not 
cross the Y-axis. Hence, there are no interior zeros for the 


corresponding eigenfunction $y of Ag. Figure 7. 
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6 (1}-o0-nx 


_ 400g 


8 {1}-° 


Figure 7. Interior Zeros of 7) 
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Recall that for each A, ; the corresponding 


eigenfunction is determined by the formula 


® .=rcos® 


Now, as A continues to increase, the _ solution 
trajectory of 6, continues to rotate in the clockwise 
direction. Eventually, the trajectory of 6, again crosses the 
line 8, (1) = o-nx when n= 1, or 6 (1) = o-x. We call this 


eigenvalue A,. Therefore, when A= 4A,, the solution trajectory 
l ] 


of 8, crosses the F¥-axis one _ time. Therefore, the 
eigenfunction & of A, has one interior zero. This process 
continues as A increases to ©, Similarly, for A<0, the 


solution trajectory of 8, rotates counterclockwise and this 
process continues as A decreases to -o. 

In summary, as A goes to +, each time that the 
solution trajectory of 8, crosses the line @(1) = otnx, the 
corresponding value for A becomes an eigenvalue of the system 
(2.28). Hence, there is an infinite sequence of A, for n= 0, 
tl, +2,... for which the second endpoint condition is 
satisfied. Therefore, there are eigenvalues, A, such that 
00 SA SAQSALS 2. and the associated eigenfunction belonging 
to A, has exactly n interior zeros in the interval -1<x<l. 


Moreover, it follows that 


limA, =© 


noe n--©& 
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In addition, any two solutions which satisfy the same 
initial conditions are linearly dependent. Therefore, each 
eigenvalue uniquely determines an eigenfunction up to a 
constant factor. This completes the proof of Theorem 5. 

For the system of differential equations (2.14), the 
discriminant of the associated Prtifer equation is positive. 


Therefore, all of the previous theorems apply. 


D. ORTHOGONALITY AND COMPLETENESS 
The objective of this section is to show that in a Hilbert 
space the eigenvectors of (2.14) are complete and orthogonal 
with respect to a weight function. [Ref. 3:pp. 344-353] 
1. Orthogonality. 
The system of differential equations (2.14) is 


equivalent to the following form 


/ 
aw) -9(2F 83) (8) «(8 3x) (F) 9 


For convenience, represent (2.34) with the notation 


(1-x?) Z,=AAZ, +BZ, (2.35) 


where 2, is the eigenvector (% ¥) associated with the 


eigenvalue 4A. A and B are the 2 x 2 matrices given below 
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~ (2x 1+3x? . (0 1+3x? 
A (2 3) and B (¢ | 


In addition, we let L stand for the linear 
differential operator, such that, when L operates on Z then LZ 


is defined by 


LZ= (1-x*) Z/-AAZ-BZ=0 


The general strategy is to find the adjoint operator, 
L,, such that Ww (LZ) - (Liw'yz is an exact derivative. Then, 
we integrate both sides to establish orthogonality. Finally, 
matrix algebra can be used to find the corresponding weight 
function. 


Let W be a 2 x 1 column vector such that 


WT(LZ) = (1-x?) W7Z’/-AWTAZ-W"BZ 
= [ (1-x?) W7Z)/- [(1-x?) W7]) ‘Z-AWTAZ-W7BZ tj 


Then, by requiring the exact derivative 


W™ (LZ) -(L°w*) Za ((1-x2) WZ], (2.36) 


it follows that 


(L°w7) Z=-{[ (1-x*) WT)! Z+AWTAZ+W7BZ) 


34 


Hence, 


L*w=-{[ (1-x”) W] +A ATW+B7W=0 


For general yp, the adjoint differential equation is 


[(1-x?) W7] ‘= -pW7A - W7B Or (2.37) 
[(1-x?) W]/ = -pAT™W - Bw 


Next, we return to (2.36). Substituting Wi (LZ) and 


(Liw')z into the equation and simplifying yields 


w? (LZ) -(L*w?) Z=—S [w?(1-x?) Z) 


which simplifies to 


wl (1—x2) Z+ (Ww? (1-x?) }/ Z=-< [w?(1-x2) Z] 


Substituting in (2.35) and (2.37) results in 


J tage. 4 pwr 4— 
(A-p) W7AZ ae le (1-x*) Z]) 


Now, we integrate both sides. This leads to the following 


orthogonality relationship 


(A-p) f° WEAZ=(W? (1-2?) Z) fi =0 Ze) 
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If the integrand is well defined, then the 
eigenvectors W are orthogonal to those of Z, with respect to 


the matrix weight function A. Next, we determine the nature 


of W. 
Expanding (2.37) leads to 
(1-x?) W=2xW-pATW-B'W (2.39) 
We write W in the form W = JE where £ is a2 x 1 colum 


vector, and J is a 2 x 2 matrix of the form 


Then, (2.39) becomes 


(1-x?) JE/=2xJZ-pA7JE -B7IZ 


Then we, premultiply both sides by zy). Note that 


giz -3. Thus, 


(1-x?) D/=2xd +pJATVJE +JB™VEZ . (2.40) 


Since 


and 
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gt | cx ox tL 
JBI (75 0 


then (2.40) simplifies to 


(1-x*) D/=pAd + (B+4xd) Z (2.41) 


where I is the identity matrix. Next, letting 


= (1-x*) -7Q and 
D=4x(1-x*) 3Q+(1-x*) -*Q 


where 2 is a 2 x 1 column vector and substituting into (2.41) 


gives 


(1-x?)Q/=pAQ+BQ (2.42) 


So, Q2 = Z where zy is the eigenvector corresponding to the 


eigenvalue wp. Back substituting results in 


W= JE = J(1-x*)7Q = J(1-x"7) "2, 
If A * wp, then the orthogonality relationship (2.38) 
becomes 


yf? gt _JSTA z 
(A-w) fo ans OO 


Since A #yp, then A-p # 0 which implies 
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f Vis JA Z, ax = 0 ° 


When A = un, then 


. JTA 
2) ee 6 (2.43) 
j, 3 (1-22)? 
Since the integrand is always positive. In addition, since 
the integrand is a polynomial, then it is finite. This 


implies Zy is square-integrable with respect to the matrix 
weight function 


JTA 
(1-x?) 2 


Initially, the (1-x?)7? term may seem to imply that the 
integrand is not well defined at +1. However, (1-x*)7 cancels 
in (2.43), and all the remaining entries are polynomials. 

Next, we determine the sequence of eigenvectors. 
Clearly, A =1, 2, 3,... all yield solutions. This is based 
on the definition of 4 such that A = n-1 where n = 2, 3, 


4,.... Now, we consider A = -l, -2, -3.... Since 


for the Legendre polynomials, the following relationships 


hold. 
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Pyar -((aea)oa) ™ *-a-2 (2.44) 


Pasa -({ (442) +1) P43 (2 . 45) 


Substituting (2.44) and (2.45) into the original 


definition of eo yields 


PP, P, P, 


a te 
n,1 n,1 
Pra Pasa P_i2 Py-3 


for A =>1. Equivalently, for A <« -l, 








This form provides the solution to the original 
differential equation (2.8). Therefore, there exists a 


sequence of eigenvectors 
© g w@ ‘ € a ¢ 
coorl ® reeerl ® ‘ly ‘ly ‘ly ‘ly reserl ®& eee 
~e -2 -1 0 1 2 » 


where (% ¥) has polynomial entries of degree < mt3. 
2. Completeness 
Finally, we determine that the sequence of continuous 


eigenvectors (8 Y), for m = 0, +1, +2, +#3,..., which are 
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orthogonal and square-integrable with respect to a matrix 
weight function, is complete. That is, we show that every 
square-integrable function F = (1-x") (f g)! can be expanded 
into an infinite series such that 


(1-24) (5) ‘> ge! (). , 


M=-0o 


Recall that there is a factor (1-x*) which cancels from both 
Sides of the equation. 
Consider the column vector (p q)' where p and gq are 


polynomials of degree < m+3. There exist constants c.. and Cy 


RB 
such that 


(a) - [ <= (2) % (¥) 


is a column vector with polynomial entries of degree <mt2. 


xats 


This follows by equating the coefficients of and solving 


the resulting system of equations for the coefficients Cs and 


c,- Since (% ¥)., and (& ¥), are linearly independent, then c., 
and c, exist. By mathematical induction, one can find 
constants ¢c, for k = -m, ..., “1, 0, 1, ... , m such that 


(5)=¢-. (#) t...4C4 (#) + (¥) +c, (¥) +--+ (#). 


Finally, we let F = (l-x"*) (€ g)! where f and g are 
any continuous functions on a finite closed interval -l<x<l, 
and let E,, €& > 0. Then, by the Weierstrass Approximation 


Theorem, there exist polynomials such that 
|(a-x2) £-P| <e, and 


|(a-x2) g-a| < @ 


for all x on [-1,1]. Hence, 


a a #2 
pm) (5) (2) (1-x*) g-qg 











< (e!] < max {e,,e,} (i) 





( (3022) 3-2) 


which implies 


(1-x*) (3) - = C, () < max {e,,e,} (7) 


Therefore, F is approximated uniformly, within an 
arbitrarily small distance € > 0, by a linear combination of 
the ($ ¥). Hence, the sequence {(S¥),} for m= ..., -i, 0, 
1, ... is a uniformly convergent sequence, which implies the 
sequence iS mean square convergent on -l<x<l with respect to 
an integrable positive weight function. Recall that (1-x*)7? 
cancels in the integrand, leaving a positive definite matrix. 
Therefore, the sequence is complete. [Ref. 2:p. 353] 
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In addition, the constants Cc, are determined in the 


following way 
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III. RESULTS AND CONCLUSIONS 


A. RESULTS 
We have shown in this thesis that there is a special 


function which takes the form 


P,(x) P,.,(x)|  (m>k) -k=0,1,2,... 


is) 
P,(X) Pas, (X) 


® (x) = 








(-1<x<1) M=Y 2,3, eee ’ 


where the elements'7 of the determinant are Legendre 
polynomials, and possesses the essential Sturm-Liouville 
properties. In addition, several interesting results are 
obtained and listed below. 
l. Hypergeometric Functions 

The Legendre differential equation is a special case 
of the hypergeometric differential equation. This suggests 
that a Beeaoiinant of the form (3.1), with other 
hypergeometric functions as elements, might possess similar 
properties. 


Recall that the recursion formula 


Paoy (x) =xP, (x) - AOD) po (x) 


and the Legendre differential equation are used extensively in 
the initial development of (3.1). In fact, it can be shown 
that since the hypergeometric function can be expressed ina 
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recursion formula of the form 


F(r+1 ; x) =A(r; x)F(lr; x) + Blr; x) P(r; x), 


then the determinant (3.1), with other hypergeometric 
functions as elements and appropriate boundary conditions, 
possesses the same Sturm-Liouville properties. 
2. Turan's Inequality 
The proof of Turan'’s inequality is a straight forward 
application of the general case (3.1). Setting n = k+l and 
substituting back into (3.1) results in a determinant that is 
of Turan's form. Recall, fromthe previous analysis of (3.1), 
that there exists a A such that (3.1) has no zeros in the 
interval -1l<x<l. Hence, Turan's inequality is established. 
3. General Systems of Differential Equations 
In general, boundary value problems of systems of 
differential equations are not well-posed. Until now, very 
little has been established concerning systems with polynomial 
coefficients. In Chapter II, several theorems are developed 
and proved for such systems. Hence, if the conditions of the 
theorems are satisfied, then the system is said to be of 


Sturm-Liouville form. 
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B. CONCLUSION 

The analysis of this one problem presents several possible 
avenues of further research. Only the first step was taken in 
a direction which could result in establishing several new and 
substantial results. For example, the next step would be to 
establish the case when (3.1) has other hypergeometric 
functions as elements. Then, one could continue the research 
on Sturm-Liouville systems of differential equations by 
examining the possibility of expanding the necessary and 
sufficient conditions for such problems. Research on this 


subject is ongoing. 


45 


APPENDIX 
I. PROOF FOR GENERAL CASE - $46) 


A. PURPOSE 

In this appendix we briefly outline the procedure for 
establishing the Sturm-Liouville properties for the general 
case, $y (). The basic procedure is very similar to that 
illustrated in Chapter II. Thus, only the important steps are 


highlighted. 


B. DEVELOP DIFFERENTIAL EQUATION 
By following the procedure in Chapter II, it is easy to 


establish the following relations for #))(x), eee and 





BY CX). 
=- (1-x?) PrPo _ PoP ([ 
gars k+1 
@/= (7k) (1=x*) pip! 4 (n-k) PLP, (2) 


(m+1) (k+1) 


PLP,+ — Ak=ntl) p pes 


oe (n-K+1) 
Y= (a-k) | AB Poa aaa CRAY Ket 


nt 


Equations (1), (2), and (3) lead to the differential 
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equation 


2 (ae) Me / 
(1 xt) | Pes Wee i (Ph? (1-x*) 3 ae a7 PxPat (ke a (Fh) |e ea 








(n-k) | (n-K+1) P,?+ oa yt mee el 


(k+1)# 


where the general solution is 


P, Pyros P, Pyot 
on ae! 
P,, Poet Or, Oni 


Applying appropriately translated boundary conditions 
eliminates all logarithms, QO: in the general solution. Now, 


setting A = n-k-1 and letting 





O-aayE (PERE and 5) 
_ 2Pr 
Ort | ere P,| - 


then (4) takes the form 





P.¢ 
2.P, (Py+ 75 Py) 


(1x2) @/- (1-2) 2 o/s (a+2) A+ = 


|@-0 5 OO 
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C. REDUCTION TO A SYSTEM OF FIRST ORDER DIFFERENTIAL 
EQUATIONS 


Next, we reduce (7) to a system of first order 


differential equations of the form 


(1-x®) @/ = ad + OF -1<x<1 
(1-x®) P/ = £O + 7? 


(8) 


where a, 8B, ¥, and Q are polynomials in x. 

First, note that Q and Q’ have no zeros in common. This 
implies that all polynomial coefficients of @, $’, and 8” are 
reduced to lowest terms and only simple zeros occur in the 
denominator. The proof is as follows. 


If Q’=0 , then either P,’ =0 or 


aap Se , (9) 





If Pi’ =0, then Q=0 only when P,=0. However, it is easy to show 
that P, and P,’ are relatively prime. Hence, when P,’=0, then 
Qe 0. 

If (9) holds, then by solving for P, and substituting back 
into Q it follows that Q=0 only when P,’=0. But, by the first 
case, this can not occur. Therefore, Q and Q’ have no zeros 
in common. 

By transforming (8) into a second order differential 
equation and equating coefficients of 8” and 8’ it is easy to 


establish that 
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y= -2x-a . (10) 


Equating the coefficients of @ and substituting (10) for y, 


results in the following equation. 


a/+ 28 __a(-2x-a) _ ad’ _ 
1-x* 1 —x* (2) 


(11) 








A+ 
~4(A+1) - = 2P,(Pyt Ph | 


Now, we isolate and examine the components of the equation 
that have Q in the denominator. The idea is to select a such 


that 





ao’ (A+1) [224 (Pet mare = 0 mod(Q) 


Substituting (6) for Q’ implies 











2 Xx Xx 2 
a| gar? ez Pht Ps) |- ara) [22x (Pt iP) | = 0 mod (Q) 


However, 
x 
Feit * Ps ore 
and © are relatively prime. This is shown by the following 
argument. 
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Set (12) equal to zero and then solve for Py. By 
substituting back into Q it follows that (12) and Q have zeros 
in common when P, and te have zeros in common. However, P, 
and Py are relatively prime. Hence, (12) and Q are 
relatively prime. 

Therefore, since (12) and Q are relatively prime, then 


there exists a polynomial M such that 


a 7 « (1-x®) 
race (A+1) P, (kK+1)2 (Py) 2+PZ ‘ (13) 
By examining (13) it follows that 


(A+1) P,+MP,2 = 0 mod (Py) 


Moreover, since P, and Be are relatively prime, then there 


exists a polynomial N such that 


(A+1) + MP, = NP, . (14) 


However, by the Euclidean Algorithm for Polynomials, there 


exist unique polynomials M, and N, such that 


M,P, + N,Py = 1 


From (14), this implies that 


M=-(A+1)M, and N= (A+1)N,. 
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Hence, there are unique polynomials M and N that determine 
a such that all the zeros of Q divide out. 


Finally, choose 8 to satisfy (11). 


STURM-LIOUVILLE PROPERTIES 


D. 
Although it is somewhat tedious, it can be shown that 
46(-Q) -[y-a]2>0 
Therefore, all of the theorems previously established 
apply. In addition, the properties of orthogonality and 


completeness follow in the same manner. 


Syl 
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